Evidence for Adiabatic Magnetization of cold Dvn Clusters 
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Magnetic properties of Dvn clusters in a molecular beam generated with a liquid helium cooled 
nozzle are investigated by Stern-Gerlach experiments. The cluster magnetizations \x z are measured 
as a function of magnetic field [B = — f .6T) and cluster size (f 7 < N < 55). The most important 
observation is the saturation of the magnetization /-L Z (B) at large field strengths. The magnetization 
approaches saturation following the power law \fj, z — fj,o\ ~ 1/vB, where /no denotes the magnetic 
moment. This gives evidence for adiabatic magnetization. 
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The first Stern-Gerlach experiments jlj on metal clus- 
ters have been performed on FeN with N=15-650 Sur- 
prisingly, in these experiments the FeN clusters were de- 
flected only in the direction of the increasing magnetic 
field, in contrast to the well known deflection patterns 
of atoms and small molecules |3|,|]]. Since then the ques- 
tion of the magnetization process of small isolated clus- 
ters has been widely discussed in literature. Khanna et 
al. suggested that the spins of the metal atoms in the 
clusters are coupled superparamagnetically resulting in 
one large spin ||. If the thermal energy is substan- 
tially larger than the magnetic anisotropy energy, the 
classical Langevin model can be applied to explain the 
magnetization of the clusters as a function of field B, 
cluster size N and internal cluster temperature T c . Ac- 
cording to this model the freely fluctuating spin relaxes 
under isothermal conditions, since the cluster vibrations 
seem to provide a heat bath which is sufficiently large 
to maintain a constant temperature during the magne- 
tization process, i.e. when the clusters enter the mag- 
netic field. Most of the experiments performed to in- 
vestigate the magnetic properties of small metal parti- 
cles were carried out at high temperatures. The results 
can be understood within an isothermal magnetization 
model using [i z = L{[1qB / '^T,.)^, where ks denotes 
the Boltzmann constant and C the Langevin function. 
Assuming T c « T„, with T n denoting the nozzle tem- 
perature, the approximation ji z = mq(mo-B /3fcsT c ) in 
the low field limit (/iqB <C fcsT c ) was applied to fit the 
magnetization curves. The obtained values for the mag- 
netic moments per atom no/N of large Con @, FeN @] 
and NIn |§ clusters with N> 2500 were consistent with 
macroscopic magnetic properties of Fe, Co and Ni. How- 
ever, the saturation of the magnetization predicted for 
the strong field limit (^oB ^> fcgT c ) of the Langevin 
model ju, z = (1 — kT c J hqB)hq which allows an indepen- 
dent determination of cluster temperature and magnetic 
moment was never observed. 

Although in many cases the Langevin model was ap- 
plied very successfully, it was found that under cer- 
tain conditions it does not explain the experimental re- 
sults correctly (7). There are two scenarios, where the 
Langevin model is likely to fail: Firstly, when the clus- 



ter temperature is so small that only very few internal 
degrees of freedom are thermally excited and therefore 
the heat bath provided for the isothermal relaxation pro- 
cess is not sufficient, and secondly, when metal clusters 
are investigated with an anisotropy energy larger than 
the thermal energy. In this case the spin cannot rotate 
freely, but is fixed to one of the crystal axis (locked spin 
model) . 

Bertsch at al. proposed an adiabatic model for the 
magnetization of clusters || with a magnetic anisotropy 
energy larger than the thermal energy of the clusters 
(locked spin model). It is assumed that no vibrational 
levels except the ground state are populated. Spin and 
rotation of the clusters are treated classically. Using sta- 
tistical mechanics, the authors derived the magnetization 
of an ensemble of classical symmetric rotors with the spin 
fixed to one of the figure axis. In the strong field regime 
(ksTji <C (J-qB) with Tn being the rotational tempera- 
ture of the cluster, the expression 
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was obtained, while in the weak field regime the equation 
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holds. At large field strengths saturation of the mag- 
netization is predicted (Equation |J). In this regime the 
independent determination of the rotational temperature 
and of the magnetic moment of the clusters is possible. 

Up to now only one experiment, performed by Dou- 
glass et al. Jl(J, has been discussed in the frame of the 
adiabatic locked spin model. The shape of the Stern- 
Gerlach deflection patterns of GdN particles with a bulk 
anisotropy energy of 10 6 erg/cm 3 at OK |Tf| generated at 
a nozzle temperature of T n = I05K were explained in 
terms of this model, but no experiments on the magneti- 
zation process have been performed so far. In this paper 
we want to investigate the magnetization of clusters far 
away from the Langevin conditions and discuss the prob- 
lem of isothermal versus adiabatic magnetization. For 
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this purpose we study cold DyN clusters (T„ = 13K) 
with a large bulk anisotropy energy (10 8 erg/cm 3 at OK 

0)- 

For the generation of DyN clusters (17 < N < 55) we 
used a pulsed laser evaporation cluster source incor- 
porated in a Stern-Gerlach molecular beam apparatus. 
The experimental setup is described elsewhere jl2|]. The 
source has been modified to produce clusters with very 
small temperatures. Now it is possible to cool the nozzle 
down to temperatures of T n =13K using liquid He. Ad- 
ditionally, the source is constructed such that the dwell 
time of the clusters in the cold nozzle channel should be 
sufficient to establish thermal equilibrium between clus- 
ters, He and nozzle. Therefore it can be expected that 
the cluster temperature before the adiabatic expansion 
equals the nozzle temperature. To favor a strong adi- 
abatic expansion through the nozzle into the high vac- 
uum of the flight tubes, which leads to further cooling 
of the clusters, a large He pressure of about lOOmbar in 
the nozzle channel is applied. In fact the velocity of the 
clusters measured behind the nozzle confirms the exis- 
tence of strong adiabatic expansion and thermal equi- 
librium between nozzle and clusters: Using the equa- 
tion Vhc = \/5RTHe/mHe flj| , where m# e denotes the 
molecular weight, R the gas constant and T# e the He 
temperature, we can calculate the He terminal speed 
Vjje, if we replace T# e by T n assuming thermal equi- 
librium between nozzle and He. At small nozzle temper- 
atures (13K < T n < 40K) the calculated He terminal 
speed Vue agrees very well with the measured velocity 
of the clusters Vc, while for higher temperatures a grow- 
ing velocity slip is observed. For example, at T!„=13K 
we find Vc — 396m/s±30m/s and Vn e — 368m/s w Vc 
versus Vc = 1200m/s±40 m/s and Vh& = 1765m/s at 
T n =300K. The fact that no slip between the calculated 
He velocity and the velocity of the clusters is observed 
for T n < 40K indicates a strong adiabatic expansion and 
a thermal equilibrium between nozzle and clusters [fl3|| . 
After being collimated the cluster beam passes the Stern- 
Gerlach magnet. The deflection is detected size selec- 
tively by a time of flight mass spectrometer in combina- 
tion with an ionization laser beam. The magnet and the 
detection unit are described in |T^| . 

For DyN clusters generated at nozzle temperatures 
T„=13K and T n =18K, we observe a shift of the Stern- 
Gerlach profile in the direction of the increasing mag- 
netic field. This indicates that relaxation processes are 
involved. Studying the field dependence of the magneti- 
zation, we obtain magnetization curves which show sat- 
uration at large field strengths as displayed in Figure |]. 

To make sure that the observed saturation of the mag- 
netization is not an effect of the relaxation time scales in- 
volved, we repeat the deflection experiments on DyN clus- 
ters with a Stern-Gerlach magnet of half of the length. 
We obtained half of the deflection and therefore the same 
magnetization. This shows that the measurement of the 
magnetization takes place under stationary conditions. 
Hence, within our resolution, there are no relaxation pro- 



cesses involved with relaxation times in the order of mag- 
nitude of the experimental time scale of typically 200/is 
which is the time DyN clusters need to traverse the Stern- 
Gerlach magnet. 

Now we turn to the question whether we expect 
isothermal or adiabatic magnetization under the exper- 
imental conditions described above. To understand the 
nature of the relaxation process it is important to con- 
sider which degrees of freedom are accessible in the clus- 
ter. As pointed out above, the clusters are in thermal 
equilibrium with the nozzle before the adiabatic expan- 
sion. After the adiabatic expansion the clusters exhibit a 
vibrational temperature T v n, which is close to the nozzle 
temperature (T v n, w T n < 18K). However, the rotational 
temperature Tr is smaller, because the adiabatic cooling 
is more effective for rotations than for vibrations [ p^[ , 
since the coupling of the rotations on the translational 
modes is stronger. To estimate the number of vibrational 
states, which are thermally accessible in DyN clusters, we 
calculate the vibrational partition sum and the occupa- 
tion number of the vibrational ground state of the dimer 
Dy 2 . As an eigen frequency we use the Einstein frequency 
of bulk Dy. We find that the ground state is occupied 
with a probability of 99.9% at T vib = 18K. Although 
it is very likely that in DyN clusters vibrational states 
with smaller frequencies than the Einstein frequency are 
available, this example demonstrates that only very few 
of the vibrational levels will be thermally accessible. To 
estimate the number of thermally accessible rotational 
states, we calculate the rotational occupation numbers 
of Dy2o ■ We approximate the complicated cluster struc- 
ture by a sphere with the density of bulk Dy and the 
mass of Dy2o- For Tr — IK we obtain an occupation 
maximum at the rotational level Jr =27 (2.2 %), where 
Jr denotes the rotational quantum number. Taking into 
account that large quantum numbers like Jr—80 are still 
populated (0.15%), the rotation of the clusters can be 
described in good approximation by the model of the 
classical rotor. 

In the next step we want to address the problem, 
whether isothermal magnetization is possible, when the 
contribution of the vibrations to the heat bath neces- 
sary for the isothermal Langevin process is neglegible, 
i.e. when the heat bath consists exclusively of rotata- 
tional degrees of freedom. This question can be answered 
by considering the entropy transfer between the spin and 
the rotational system employing the Second Law of Ther- 
modynamics. According to Boltzmann the entropy S 
can be written as S = — fc^Wi In u>i, where cj, de- 
notes the occupation probability of the state i Q . The 
entropy loss ASj provoked by the magnetization pro- 
cess is due to fixing the orientation of the total angular 
momentum J of the clusters in the magnetic field di- 
rection z at large field strengths. This corresponds to 
the saturation of the magnetization. Without field the 
cluster has 2 J + 1 possibilities of equal probability iOi 
to orientate the total angular momentum versus the z- 
axis. Since we observe experimentally the saturation of 
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the magnetization (see Fig. |l|), the order of magnitude 
of the magnetic moment /j,q = gj^J J(J + 1)/xb> where 
gj denotes the g- factor of the cluster, can be estimated 
by taking into account that the saturation magnetization 
(is is approximately equal to the magnetic moment 
Assuming that the g-factor of the clusters is in the or- 
der of magnitude of the g-factor of the ground state of 
the Dy atom gj — 1.33 ~ g,j, one obtains for Dy2o with 
fi s « /iQ rj 9/i B a total angular momentum of J rj 6 
and hence an entropy loss ASj = J2i=i hi(l/13)fcs/13 = 
— fcslnl3. According to the Second Law of Thermody- 
namics the rotational entropy ASr > —ASj has to aug- 
ment at least by the same amount. This entropy gain 
changes the rotational temperature Tr according to the 
equation ASr = cvrATr/Tq, where cvr is the heat ca- 
pacity and To denotes the temperature of the clusters, 
before the field was applied. We have shown above that 
many rotational levels with large quantum numbers are 
occupied in Dy2o ■ Therefore we approximate cy r by the 
high temperature value cvr = 3fcs/2. Hence the temper- 
ature change due to the magnetization process is given 
by the expression AT R /T > ln(2J + 1)2/3. For clus- 
ters with an initial rotational temperature To « IK the 
temperature rises during the magnetization process by 
ATr w 1.7K. This significant change of the rotational 
temperature shows, that in our case it is not suitable to 
use the isothermal magnetization model. 
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FIG. 1. The magnetization curve of Dy2o is shown. The 
solid line is the fitting function according to the adiabatic 
model by Bertsch et al. in the strong field limit [9]. The 
dotted line is the magnetization curve expected in the weak 
field limit 

Now let us turn to the question, whether adiabatic 
magnetization is possible. For the magnetization pro- 
cess the time t s « 10~ 5 s needed by the clusters to move 
from a zero field region into the magnet with the full field 
applied, is important, because adiabatic magnetization is 
only possible when the magnetization takes place at much 
larger time scales than the relaxation process. Since the 
typical time scale for cluster rotations is 10 -9 s, the mag- 
netization can take place adiabatically. Hence, we con- 
sider whether the adiabatic model proposed by Bertsch 
et al. can be applied to evaluate our experimental data. 



The assumptions made in this model (classical rotation, 
classical spin and locked spin) reflect our experimental 
conditions. It has been discussed above that the rotations 
can be treated classically. In first approximation the spin 
of the cluster can be treated classically as well, because 
the saturation magnetization of the clusters indicates the 
existence of large total angular momenta (see Dy2o). The 
locked spin model is justified, since the thermal energy 
of the clusters (T c < 18K) is much smaller than the mag- 
netic anisotropy energy. Using the anisotropy energy of 
bulk rare earths (« 10 8 erg/cm 3 at OK Q), we estimate 
for Dy2o an anisotropy energy which corresponds to a 
temperature of about 500K. Therefore we use the adia- 
batic model proposed by Bertsch et al. to fit the mag- 
netization curves of DyN clusters, recorded for N=17-29 
at T„=18K and for N=32-55 at T„=13K. Since most of 
the data measured at T„ = 18K and at T„ = 13K belong 
to the strong field regime (see Fig. |l|), we use Equa- 
tion |l| to fit our magnetization curves. After plotting 
the magnetizations versus 1/VB, the magnetic moments 
and the rotational temperatures Tr of the Dy^ clusters 
are determined by linear regression, as shown in Figj^. 
In average Equation [j] gives a good fit for the clusters 
generated at T„ = 18K and at T„ = 13K. Only mag- 
netizations measured at very small field strengths, i. e. 
5=0. 22T for T n = 18K do not fit the expected behavior, 
because they do not belong to the strong field regime. Tr 
is calculated from the slope of the line fit and depends 
on the prefactors imposed by the structural assumptions 
and the assumptions about the distribution of the rota- 
tional energy. Hence the temperature scale depends on 
the geometrical details of the clusters, while the magnetic 
moment is independent. 
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FIG. 2. Magnetizations of Dy2o plotted versus B~ 1/2 . The 
solid line denotes the fitting function obtained by linear re- 
gression according to the adiabatic model. 

Figure [5]a shows the magnetic moments per Dy atom 
Ho/N of DyN clusters with A=17-55 and their rota- 
tional temperatures Tr obtained by applying the adia- 
batic model as described above. In average the magnetic 
moments per atom range between 0.3 to 0.6 fiR- In com- 
parison to the magnetic moment of Dy \id v = 10.6/xb in 
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the ferromagnetic bulk phase the magnetic moments of 
the clusters are smaller by a factor of 20. This suggests 
that the magnetic ordering in the clusters is rather anti- 
ferromagnetic than ferromagnetic (J = 52* J» ~ 0, with 
J being the total angular momentum of the cluster and j 
being the total angular momentum of the Dy atom), al- 
though the rotational cluster temperatures are well below 
the Curie temperature of the bulk Tc(Dy)=86K. Keep- 
ing in mind that the ferromagnetic ordering in rare earth 
metals is determined by indirect coupling through the va- 
lence electrons (RKKY interaction [|l5| ) , this result is not 
very surprising and similar behaviour has been observed 
for other rare earth clusters JL6| |. Since the structure of 
small clusters differs strongly from the bulk lattice to 
compensate surface effects, the wave functions of non lo- 
calized electrons like valence electrons change in response 
to the change of the long range structural ordering. As a 
result the coupling between the 4f 9 cores of the Dy atoms 
in the clusters differs from the coupling in the bulk, since 
electrons in irregularly shaped cluster orbitals are polar- 
ized instead of electrons in orbitals described by regu- 
larly oscillating Bloch functions. The theoretical study 
performed by Pappas et al. on Gdi3 illustrates the ef- 
fect of the modified exchange coupling on the magnetic 
structure of small rare earth particles very well In 
Figure ||b the rotational temperatures Tr of the clusters 
are shown. Since the rotational temperatures (0.2-2. OK) 
are much smaller than the nozzle temperatures (T„=13K 
and 18K), we conclude that the adiabatic cooling of the 
rotational degrees of freedom after the expansion takes 
place very efficiently, as it is suggested in Q. 

Although the adiabatic model developed by Bertsch 
et al. fits the recorded magnetizations very well in the 
strong field limit, the overall shape of the magnetization 
curve does not match the theoretical prediction in Ref. 
H for a spherical cluster. The theoretically predicted 
adiabatic magnetization curve shows a linear dependence 
of the magnetization at small fields (Equation |2[) like the 
isothermal magnetization, whereas the experimental data 
(Fig.|l|) suggest the existence of higher order terms. Since 
it cannot be assumed that all DyN clusters in a size range 
of N= 16-55 are spherical, the breaking of the spherical 
symmetry could probably account for the differing curve 
shapes at small fields, as it was shown by quantumme- 
chanical calculations [18|. 

In summary, the magnetization of DyN clusters gen- 
erated at low nozzle temperatures can be understood by 
assuming adiabatic magnetization, while an isothermal 
process is not consistent with the experiment. By mea- 
suring the saturation of the magnetization, the magnetic 
moment and the temperature of the clusters can be de- 
termined at the same time independently. The rotational 
temperatures obtained suggest a strong adiabatic cool- 
ing. The DyN clusters exhibit a nearly antiferromagnetic 
spin order, although bulk Dy is ferromagnetic in the tem- 
perature range studied. This can be understood in the 
frame of the RKKY theory. The complicated size de- 
pendence of temperature and magnetic moment of the 



clusters remains an open question. 
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FIG. 3. In (a) the magnetic moments per Dy atom /xo /N of 
DyN clusters are plotted versus N, in (b) the rotational tem- 
peratures Tr. Filled squares denote DyN clusters generated 
at T n =18K, empty squares denote DyN clusters generated at 
T„ = 13K. 
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